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|/-\ . We consider circular and elliptic quantum dots with parabolic external confinement, containing 

f*"^ ' — 22 electrons and with values of r^ in the range < r^ < 3. We perform restricted and unrestricted 

f~^ ^ Hartree-Fock calculations, and further take into account electron correlations using second-order 

fSJ ■ perturbation theory. We demonstrate that in many cases correlations qualitatively change the spin 

structure of the ground state from that obtained under Hartree-Fock and spin-density-functional 
j^ I calculations. In some cases the correlation effects destroy Hund's rule. We also demonstrate that the 

correlations destroy static spin-density waves observed in Hartree-Fock and spin-density-functional 

calculations. 
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I. INTRODUCTION 



(-H ' Recently many theoretical and experimental studies have focused on quantum dots, in which electrons are confined 

I ' to a mesoscopically small region. These structures are of interest both for studying fundamental physics and for 

Q\ ' potential nanotechnology applications. With advances in materials science, semiconductor quantum dots have become 

particularly popular. Most semiconductor dots consist of layers of semiconductor with the electrons confined by the 

band structure to a layer interface, forming a quasi-2D electron gas. The electrons are further confined in the plane 

to form the dot, usually by charged gates. Quantum dots containing many electrons are often referred to as "large 

dots" , and those with few (often N < 20) as "small dots" . A large dot has a relatively dense electron spectrum and 

therefore poorly controlled external parameters like the shape of gates, impurities, etc. bring an intrinsic external 

i-rt " random component into any description of the dot. For a small dot the uncertainty in the parameters is negligible. 

P5 The external confining potential can be approximated as Uext = hiTT'*{^x^'^ + ^yy'^)i where m* is the effective mass 

O of the electron in the semiconductor, and for circular dots uJx = ^y Such small dots are in essence artificial atoms. 

I ^i . However, they have the advantage of allowing direct measurement of transport properties, and experimental control 

over basic parameters such as the size and strength of the confining potential and the number of electrons in the 

^^ dot. They can also be driven by external magnetic fields that would be inaccessible in atomic physics with realistic 

k>" ' magnetic fields. 

^^ \ Small circular quantum dots have been investigated experimentally by Tarucha et al} and by Kouwenhoven et 

^Ll • al.^ . These studies have clearly demonstrated the shell structure of the energy spectrum, the effect which one would 

i(»^ [ expect from analogy with atoms or nuclei. Where there is a shell structure we might also expect Hund's rule to be 

^^ . obeyed for the total spin of the dot. Quantum dots have attracted the attention of numerous theoretical studies, for 

■n^ ' reviews see Refsi^. The electronic structure of 2D small quantum dots has been studied using spin density functional 

C^ , theorj&SiSiifliiLi^ as well as by restricted and unrestricted Hartree-Fock methodsi^. These studies have also clearly 

^i ■ demonstrated the shell structure of circular dots and they have confirmed the validity of Hund's rule for such dots. 

2 ! As expected, elliptical deformation destroys both the shell structure and Hund's rule. The spin-density-functional 

^ ■ calculations also lead to spin-density wave states. However, as has been pointed out by Hirose and WingreeniS such 

'■ ' states are artifacts of broken spin symmetry in density-functional theory. This is quite similar to what has been 

Jh • known to occur in unrestricted Hartree-Fock since the work of Overhauser^^. 

Q ' As mentioned, there are similarities between small quantum dots and atoms. Therefore theoretical methods used 

Q , in atomic physics could be efficient in studies of quantum dots. It is known that for few (2-4) electron atoms the 

ILJ best available method is the pair equation method. The configuration-interaction method also works well in this case. 

. ^H ', However for multi-electron atoms the best available method is Hartree-Fock with further account taken of correlations 

J^ ■ using many-body perturbation theory. This method is used nowadays in practically all precise calculations for multi- 

H ' electron atoms, see e.g. Refa^SiSiS^, but to the best of our knowledge it has never been used before for quantum dots. 

. . . In the present work we perform such calculations. Our results for total energies and shell structure are rather similar 

to already known ones. However the results concerning the spin structure are different. We never observe a spin 

density wave in a state with total spin zero. The wave can appear at the Hartree-Fock level, but taking correlations 

into account restores the rotational symmetry. So, on this we agree with work&10»iS and disagree withSii^. Studies 

have also been done in which the rotational symmetry has been explicitly restored after unrestricted Hartree-Fock via 

projection techniques^ 

In the circular dot we have found that the total spin sometimes varies from that predicted by Hund's rule. For 
iV = 10, 15, 16, and 17, we found that the ground states have minimal spin as opposed to the previous result of 
maximal spin within the confines of Hund's rule. We have also studied elliptically deformed small dots, paying special 



attention to the case ujx/ojy = 2. It is well known from nuclear physics^^ that the shell structure is restored when the 
ratio of frequencies is a rational number. This is certainly true for a parabolic potential. For a quantum dot the total 
self-consistent potential is not parabolic (even if the external potential is parabolic), but nevertheless we show that 
there are some peculiarities in the spectrum of the dot at Lo^jujy — 2. 

The structure of the paper is as follows. In Section^ we present simple semiclassical estimates which allow us 
to relate the number of electrons, the strength of the confining potential, and the relative strength of the Coulomb 
interaction. Section IIIII is devoted to an overview of the computational method. Section IIVI concerns the two- and 
six-electron problems in the external parabolic potential. The two-electron problem can easily be solved exactly, 
while for the six-electron problem we compare to the configuration-interaction study of Reimann et al^. Comparing 
these accurate energies with our restricted and unrestricted Hartree-Fock results we demonstrate how unrestricted 
Hartree-Fock generates the spin density waves and how correlations restore the rotational symmetry. Based on the 
results of this section it appears that the method gives reliable results to at least rj < 3, although there appears to be 
some N-dependence in the accuracy. Finally, Section presents our results for circular parabolic dots, and Section 
IVII presents results for elliptical dots. 

II. ANALYTICAL ESTIMATES OF PARAMETERS 

In this section we use atomic units for energy and length 






IO2A « lOnm . (1) 



Here e is the elementary charge, m* is the electron effective mass, and e = keq is the dielectric constant of the material. 
Numerical values are presented for m* = 0.067me and e = 12.9eo, which correspond to the GaAs commonly used in 
quantum dots. The confining external potential is of the form 

TT -'^/22| 2 2\ -'^22 /r,\ 

Uext^^iuJ^x +ujyy)^-ujr . (2) 

Assuming that the number of electrons in the dot is large, iV 3> 1, we can use the standard semiclassical expression 
for number density in a 2D circular dot of radius i?^^ 

n{r) ^noy/l-r^/R^ , (3) 

where no = '6N /{2ttR^). The electrostatic potential energy at the origin and at the edge of the dot is 

U{R) = ^noR. (4) 

In the Thomas-Fermi approximation the number density at r = is related to the Fermi momentum at this point 

p\ = 2nno , (5) 

and the self-consistency equation reads 

At A^ ::^ 1 the Fermi energy is small (ex N~-^^^) compared to the self-consistent potential. So Eq. © implies that 
{Uext + C^)lr=i?, = {Uext + C^)lr=o- "^^^^ gives the standard semiclassical relation between the size of the dot and the 
number of electrons (see for example Refii^) 

R = rsV 
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no = ( ^ ) ujiN^ . (7) 



The number density of electrons varies from point to point. Therefore the interaction is less important at the centre 
and more important towards the edge of the dot. To characterize the interaction we define an average r^ according 
to the equation 



vrrf TV = Tri?^ , or i? = TsWN . (8) 

Then, using Q one finds 



uj-3N~s . (9) 



Note that this definition of r^ differs from that accepted in Rcf."* by a factor of (37r/4)^ k, 1.33. This difference is due 
to the semiclassical estimate we use for the dot density in Eq. 0, as opposed to the assumption^ that the density is 
roughly flat. 

The total energy of the dot scales as 

EozNir-^ . (10) 

We will use Eqs. Q, ((SJl, ®, and (|10|) not only for circular dots, but for elliptic dots as well. In this case we will 
take uj — ^ijjx^y, R = \ab. 

We stress once more that Eqs. {TJ, (0, (O, and (|10|l are written in dimensionless atomic units. Using Eq. (Q 
one can easily convert them to other units. For example a GaAs quantum dot with r^ = 1 and containing about 20 
electrons should have a radius of the order of i? « 500A « 0.1/^m. The confining potential hw must be on the order 
of 8 meV once any screening effects from the surrounding layers are considered. 

III. METHOD 
A. The model 

The electron wavefunctions are modeled on a square lattice. Essentially we consider an Anderson model with 
Hamiltonian 



H '^ 



-{Ujlxf + UjIv1)\ c\„C^a - Yl ^L^ja + 3 H 1^.. 1 ^. I ^aC^acl^Q/g ■ (H) 
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Here c]^ is the creation operator of an electron with spin projection a (t or |) at site i of the two-dimensional 
square lattice with lattice spacing equal to one. (ij) represents a sum over nearest neighbor sites. The single particle 
dispersion corresponding to the Hamiltonian (|11|) is 

Ep = 4 — 2 cos Pa- — 2 cospy ~ p^ . (12) 

At larger momenta the dispersion deviates from quadratic dispersion, so to simulate the real quadratic dispersion 
we need to have many lattice points within one period of the electron wave function, which is on the order of one 
atomic unit. For this reason, the units of length used in the Anderson Hamiltonian Hill) are not atomic units. The 
Hamiltonian is written in units in which the lattice spacing is one, where the lattice spacing needs to be significantly 
smaller than one atomic unit. 

The Coulomb interaction q"^ /{vi — rj\ is singular at i = j. This is an unphysical singularity as, firstly, if the integral 
were performed in continuous instead of discrete 2D space there would be no divergence, and secondly as in a real 
dot there is transverse (orthogonal to the plane) confinement with some finite effective length, typically on the order 
of lOnm. 

Here, we avoid adding a finite width to our model in order to be able to compare our results to exact and 
configuration-interaction 2D calculations. Instead we consider that when r.^ = r^, the effective distance between 
electrons should not be taken to be zero, but some finite value characteristic of the size of the lattice cell. We obtain 
an approximate value for the integral at r^ = r^ by treating the wavefunctions as constant over the unit cell and 
taking the integral of q^ /\ri — rj\ over the cell. That is, when Fj = Vj, we take 

drdr' w 0.34, (13) 

lattice cell I'-" ~ '■" I 



where the "inverse characteristic distance" of 0.34 is in units determined by the lattice spacing. The Hartree-Fock 
results are slightly sensitive to the value given here. However, when correlations are taken into account the results 
are not sensitive to this value. 

As a first approximation for the ground state of the Hamiltonian 111|) wc use the unrestricted and restricted Hartree- 
Fock (HF) methods, modified according to the Optimal Damping Algorithm of Cances and Le Brisi^. Recall that the 
HF equations are of the form 

-ff.T.V'j =e^V'^ , (14) 

where V'i(^) is a single electron orbital which has energy e^, and spin projection on the z-axis at — ±^. The HF 
Hamiltonian matrix is of the form 

(r|iJ,J/) = {r\Hs\r') + {r\Ud^r\r') - {r\Uexch,aAr') , (15) 

where 

{r\Ud^r\r ) = Srr' q }^ }^ 7^—^. , 

ek<eF r" ' ' 

{r\U^..H,.. W) ^q'Y. -5.... '^'f^'^f^ , (16) 

z^ r-r' 

are direct and exchange interactions. The single-particle Hamiltonian has only diagonal and nearest neighbour nonzero 
matrix elements, 

{r\Hsy) = 6rr' U + ^{uJIx' +Zjy)] - S^rr'> ■ (17) 

Hear 5rr' is the usual Kronecker delta. 6^rr'> = 1 if r and r' are nearest neighbours and (5<r-r'> = otherwise. The 
total energy of the system is 

In the unrestricted HF (UHF) method the single particle orbitals with spin up are completely independent of those 
with spin down. Therefore, generally UHF spontaneously violates the rotational invariance of the Hamiltonian l|ll(l . 
The total spin S does not commute with the UHF Hamiltonian and hence only the z-projection of the total spin, Sz, 
is conserved. To characterise this effect we introduce the parameter "spin separation" 

An=i/"|n|(r)-ni(r)|dr . (19) 

Here n^ is the total density of electrons with spin up and n^ is the total density of electrons with spin down. The 
spin separation measures the spatial separation of up and down electrons. It can be nonzero even for a state with 
Sz = 0. Strictly speaking the spontaneous violation of rotational symmetry and hence nonzero value of An for a state 
with 5z = is not a physical effect, but a byproduct of the method. Nevertheless, a calculation of An can shed light 
on the physics of the system. For example, zero An indicates that there is no spontaneous symmetry violation. If 
An « 1 in the sector with Sz — there are two possibilities: cither the total spin is S* = 1 and we effectively see the 
state |5 = 1, 5*2 = 0), or the state is |5 = 0, 5^ = 0). By performing the UHF calculation in the sector with 5^ = 1 
one can try to distinguish between these possibilities. 

In the restricted HF (RHF) method the HF equations are exactly the same, but in addition we impose the constraint 
that electron orbitals with spin up and spin down are identical. Hence RHF automatically respects the spin rotational 
invariance and describes the state with 5* = 0. The RHF total energy is always higher than the UHF total energy. 

B. Correlation corrections 

The leading correlation correction to the total energy is given by the diagrams in Fig. ^ The corresponding formula 
reads 





FIG. 1: The leading correlation correction to total energy 
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where the Coulomb matrix element is defined as 
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(m,n|y|i, j) = (5^„.^.(5^„a, ^V'm(r)V'«(rO |^_^,| V'»(r)V'3(i'0 • (21) 

r,r' 

It is well known that in a 3D uniform electron gas the direct diagram in Fig. ^is logarithmically divergent due to 
the long-range nature of the Coulomb interaction. Hence it requires consideration of screening effects (higher orders 
of perturbation theory) even at very small TgS^ Fortunately for the 2D case the diagram is convergent. Moreover we 
consider a relatively small dot and in this case the second order correlation diagram is convergent even in the 3D case. 

C. Restricted HF with an odd number of electrons 

There is no difference between an even and an odd number of electrons when performing UHF calculations. However, 
a restricted Hartree-Fock calculation with an odd number of electrons is not straightforward as the spin-up and 
spin-down electrons must feel a different exchange potential, breaking the symmetry required by RHF. The Slater 
determinant of the problem consists of {N — l)/2 doubly occupied orbitals which represent a core with total spin zero, 
and also one upper orbital which contains only one electron. To preserve the rotational invariance of the core we use 
the same method used in atomic physics; we average over polarizations of the unpaired electron. This means that we 
introduce occupation numbers rii where Ui = \ for any core orbital, ne| = riej, = 2 where e is the external orbital, 
and n^ = otherwise. The Hartree-Fock equations (|14 |l - H18f) are modified appropriately. For example instead of Eq. 
(|18|l . the total energy reads 

■E^RHF = ^ [niti - -ni{Udir)i + -zni{Uexch)i\ ■ (22) 

i 

This procedure defines our first approximation and it corresponds to the many-body state with total spin S = ^. 

The next step is to take into account the fact that the real state contains one |e t) electron and no |e |) electron. 
We use perturbation theory, which automatically guarantees that the total spin remains S — ^. The transition from 
two external "half-electrons" to one "whole" external electron results in a perturbation acting on the external electron 
AVeT ~ ~5Dir(e)j and perturbations acting on the core electrons ISVi^ — — iEx(e) and ISVii — 5Ex(e). Here Dir(e) 
and Ex(e) represent direct and exchange interaction with the "external" orbital ^e- Applying these perturbations one 
finds the following first- and second-order corrections to total energy (|22|l . 



<5i5(i) = -i(e,e|F|e,e 



4 
^^(2) ^ ly |(e,e|F|fc,e)p ^ I yy \{i,e\V\eM^ ^ I yy \{i,e\V\e,k)\\ (33) 

k>e i<e k>e i<e k>e 

We use the notation defined in Eq. H21|l . Note that (|23|l is not the same as the correlation correction. The correlation 
energy has to be calculated separately using (|20|) . In H23|) we go to second order in order to retain formal accuracy 
comparable to the correlation terms H2U|) . However, in practice we have found that SE^'^' is always very small and 
can be neglected. 





FIG. 2: Exact (left) and Hartree-Fock (right) charge densities (atomic units) for a two-electron dot at r^ — 0.77. In this 
case restricted and unrestricted Hartree-Fock results are practically identical. Note that r^ here is given by the semiclassical 
estimate made in Section^and, for only two electrons, will not be exactly equal to -Ecouiomb/£'Kinetic. 





FIG. 3: Exact (left), restricted Hartree-Fock (central), and unrestricted Hartree-Fock (right) charge densities (atomic units) 
for a two-electron dot at rs = 4.87. Note that rs here is given by the semiclassical estimate made in Section |H| and, for only 
two electrons, will not be exactly equal to -Bcouiomb/£'Kinctic. 



IV. TWO PARTICLES IN A CIRCULAR PARABOLIC POTENTIAL 



The two-body problem can be solved exactly. The Hamiltonian of the problem is 



«-|^ + li + ^'"">'i^ + i-i') 



|ri -r2 



(24) 



Using a centre-of-mass coordinate and a relative position coordinate: 

1 



R 

r 



(ri+ra) 



2 

ri -r2 



one can transform the Hamiltonian to 

H 



4:771 Oil-' m or^ 4 r 



(25) 



(26) 



Variables in the corresponding Schroedinger equation are separated and hence one can easily find the exact solution of 
the problem by solving the radial equation for r numerically. The ground state has relative orbital angular momentum 
equal to zero, I = 0. Therefore, because of Fermi statistics, the spin of the ground state is also zero, S = 0. The first 
excitation has I = S = 1. For clarity, in Fig. 0]we plot values of Erg versus r^, where E is the energy in atomic units 
and r-g is defined according to lO at A^ = 2 (see commenli^). 

Fig. ^ shows exact energies for S* = and S = 1, UHF energies for Sz — and 5^ == 1, and the RHF 5 = energy. 
Fig. ^p shows the same energies but with second order correlation corrections included H2UI) . In this case the dashed 
lines show E'uhf + i^-^uhf at 6*2 = and Sz = I- 

In Fig. [3 and Fig. |21we show exact, RHF, and UHF electron densities in the ground state for r^ = 0.77 and 
rs = 4.87 respectively. At small r^ both the restricted and unrestricted {Sz = 0) Hartree-Fock methods give results 
very close to the exact S = solution. However, at larger Vg the ground state UHF energy is lower than that obtained 
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FIG. 4: Total energy versus Vs for the two particle problem. We plot the total energy in atomic units multiplied by r^. Exact 
solutions for S — and 5=1 are shown in both Fig.a and Fig.b by solid lines. Fig. a also shows the UHF energy with 
Sz = and Sz ~ 1 (the dashed lines) and the RHF energy (the dotted line). Fig.b shows the same energies with second-order 
correlation terms 1121)1 included in all cases. Lines joining data points are given as a guide to the eye only. 

in the restricted Hartree-Fock method (see Fig. 0^). The UHF energy for the ground state approaches the energy of 
the first excitation with S* = 1 at larger values of r^ . 

This illustrates that the UHF method spontaneously violates the spin and rotational symmetry of the problem and 
mixes the true ground state, which has 5 = 0, with the state \S — l,Sz ^ 0)- At rs > 2 the 5 = 1 component 
dominates. This leads to early "Wigner crystallization" which is a byproduct of the UHF method. This is confirmed 
by the pictures of electron density at r^ = 4.87 shown in Fig. |31 There are two separate peaks in the UHF electron 
density; one corresponds to the electron with spin up and another to the electron with spin down. The UHF spin 
separation parameter H19() is An = 1.7. At the same time neither the exact solution nor the RHF solution indicate 
any spin separation. 

Taking correlations into account drastically changes the situation. The total RHF+correlation energy shown by 
the dotted line Fig. ^ is close to the exact energy up to very high r^ for the two-electron problem and is below the 
corresponding energy obtained with the UHF-|-correlations method. Thus the inclusion of correlations leads to the 
choice of the RHF solution as the ground state, hence restoring the rotational invariance at 5 = 0. Certainly at very 
large rg, all Hartree-Fock methods fail. 

Fig. [3 gives a similar comparison for the six-electron dot, with configuration-interaction values taken from the 
study of Reimann et ali^. Here Fig. [SJ^ shows CI energies for 5 = and 5 = 3, UHF energies for Sz = and Sz = 3, 
and the RHF 5 = energy, while Fig. ^ shows the same energies but with second order correlation corrections 
included 1^0)1. We see the same behaviour as in the two-electron problem, although the third and higher order terms 
not included in our method are larger when N=6. Again the true ground state is 5 = 0, and as r^ increases the UHF 
Sz = solution breaks the symmetry of the problem and converges to the polarised solution. By using RHF as a 
base for the perturbation theory we remove the spin density waves and restore the rotational invariance of the ground 
state. 

We point out that the UHF and the UHF+correlations method both work pretty well for the polarised 5=1 
and 5 = 3 states, see Fig. 0]i,b. This is because the coordinate wave function is antisymmetric with respect to the 
permutation of electrons and hence the electrons are well separated in space. 



V. THE MULTI-ELECTRON CIRCULAR QUANTUM DOT 



The second difference of the total energy is defined in the usual way 

A2(A^) = E{N + 1) + E{N - 1) - 2E{N) 



(27) 



where N is the number of electrons in the dot. Fig. |B^ shows calculated values of A2 against N in the round dot 
with confinement frequency Wat = 0.73 (atomic units) which corresponds to rg = 1.48. . .0.99 over the range of N 
(see Eq. ISJ. We also find the ground state spin. Values of the spin at even values of A^ are marked above the lines. 
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FIG. 5: Total energy versus rs for the six particle problem, with CI energies taken from the configuration-interaction study of 
Reimann et ali~ We plot the total energy in atomic units multiplied by rg. Exact solutions for 5* = and S — 3 are shown in 
both Fig.a and Fig.b by solid lines. Fig. a also shows the UHF energy with Sz = and Sz — 3 (the dashed lines) and the RHF 
energy (the dotted line). Fig.b shows the same energies with second-order correlation terms (I20II included in all cases. Lines 
joining data points are given as a guide to the eye only. 



Fig.a 



Fig.b 




3. 

< 



10 15 

N 




FIG. 6: The circular dot with confinement frequency o^at ~ 0.73 which corresponds to rs — 1.48 . . . 0.99 over the range of A'^. 
a: A2 against A'' from unrestricted HF (bottom line) and from HF-fcorrelations (top line). At even values of A'^ we also show 
the spin of the ground state. 

b: The excitation gap with spin change against A'^ (correlations are taken into account). At each N we also show the spin of 
the ground state and the spin of the excited state. 



The bottom line in Fig. |S^ describes naive HF calculations with no correlations included. These are all UHF results 
since the raw UHF energy will always be lower than the raw RHF energy. To determine total spin of the UHF state 
we perform calculations for different values of Sz and set S to the highest of the Sz which give degenerate energies. 
In cases when S ^ 0, we see very degenerate energies up to the maximal spin, followed by a much larger energy gap 
for Sz > S, and this indicates that UHF almost preserves the rotational invariance. These UHF results agree with 
those obtained previously using spin density functional theory'*'^. The 2D parabolic shell structure is reflected in the 
A2 values. We see peaks at A^ = 2, 6, 12, and 20, consistent with the shell model ls^2p'*3s^3(i'*4p^4/'*. As one would 
expect, the Hartree-Fock spin structure agrees with that predicted by Hund's rule. 

The top line (offset) in Fig. |B^ shows results obtained by taking correlations into account. At each N we perform 
UHF-|-correlation calculations for different values of Sz and we also perform a RHF+correlation calculation for 5 = 
(or S = 5 if A^ is odd). This allows us to determine the ground state energy and spin. Fig. |SJl shows that correlations 
influence the values of A2 for some N, but retain the shell structure. For some N where the shell is open, i.e. N — 10, 
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FIG. 7: The circular dot with confinement frequency Usx — 0.26 which corresponds to Vs — 2.91 . . . 1.95 over the range of A''. 
a: A2 against A'^ from unrestricted HF (bottom hne) and from HF+correlations (top line). At even values of TV we also show 
the spin of the ground state. 
b: The excitation gap against N as defined by -Euhf, Sz =1 — Eb.hf, s=o (correlations are taken into account). 



15, 16, 17, the correlations also change the spin from 1, 



I to 0, 



0, 



In these cases, including correlations 



leads to a violation of Hund's rule. However, in these cases we see that the excitation gap with change of spin 

A = Ei{N) - Eq{N) 



(28) 



is very small in these dots. The spin gap becomes large only for N corresponding to closed shells. Values of A versus 
N at Wat = 0.73, with correlations taken into account, are plotted in Fig. ^p. Unfortunately our technique does not 
allow the calculation of the gap without spin change. Note that this gap is not equal to the difference of the single 
particle Hartree-Fock energies. We should point out that at iV = 16, where the calculated excitation goes from S — 
to S = 2, we cannot accurately calculate the energy at S* = 1 as only Sz is conserved in the HF method. This means 
we cannot observe the 5* = 1 excitation - most likely it has energy comparable to that of the S = 2 excitation. 

Fig. [3 shows the same quantities A2 and A for dots with weaker confinement tJat = 0.26. This corresponds to 
stronger interaction, r^ = 2.91 . . . 1.95 over the range of A''. In this case correlations substantially influence the values 
of A2, although the shell structure on A2 is destroyed by electron-electron interactions whether or not correlations 
are included. It should be noted that in the UHF case here it becomes difficult to determine which energies should 
be treated as degenerate and which should not. The UHF spin values given on the bottom line of Fig. [7^ should 
be taken with caution. This is in contrast to Fig. |B^ where the degeneracies are clear-cut. When correlations are 
included, the ground state spin becomes minimal for almost all N. We see in Fig[7|3 that the "spin gap" of Eq. I|28|) 
grows linearly with N. At the same time, the UHF+correlations S^ = 1 calculation in this case generates a highly 
excited state with a strong spin density wave (the spin separation parameter is An ^ 4 — 5) . This implies that the 
accuracy of the UHF-Fcorrelations method drops dramatically at the larger N due to the developement of unphysical 
spin density waves, and that the gap in Fig. Et* is not a true spin gap. We believe that the S = ground state, as 
determined by RHF-|-correlations, is fairly accurate. In principle, based on this ground state, one could find the first 
physical excitation using the time-dependent HF-f-correlations method. However this is a very involved calculation 
and beyond the scope of the present work. 



VI. ELLIPTICAL DOTS 



Studies of elliptical dots are interesting from a theoretical point of view. Also, an elliptical deformation is a first way 
to model deviations from circular geometry, which always exist in real dots. A lateral dot, for instance, has contacts 
at the sides through which the electrons tunnel onto the dot, and these alter the shape of the confining potential. 
Previous numerical studies have found that a deformation of the dot destroys the shell structure. ^'®'^° In the present 
work we investigate a simple deformation of the dot by taking ujx > ujy- 

Fig. |SJi shows calculated values of A2 against N for a dot with the small ellipticity uJx/oJy = 1.1 and with 
(Jat = ^'^x'^y — 0.73, which corresponds to Ts — 1-48 . . .0.99 over the range of N . The value of the ground state spin 
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FIG. 8: The dot with eUipticity Lo^jiOy = 1.1 and with Uat = ^JuJ^UT^ — 0.73 which corresponds to r^ — 1.48 . . . 0.99. 

a: A2 against A'^ from unrestricted HF (bottom line) and from HF+correlations (top hne). At even values of A*' we also show 

the spin of the ground state. 

b: The excitation gap with spin change against A'^ (correlations are taken into account). At each N we also show the spin of 

the ground state and the spin of the excited state. 
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FIG. 9: The dot with eUipticity LOxjiOy = 2 and with tJat = y/uixUjy = 0.73 which corresponds to rs = 1.48 . . . 0.99. 

a: A2 against A'^ from unrestricted HF (bottom line) and from HF+correlations (top line). At even values of A' we also show 

the spin of the ground state. 

b: The excitation gap with spin change against A'^ (correlations are taken into account). At each N we also show the spin of 

the ground state and the spin of the excited state. 



at even values of N is marked above the lines. The bottom line in Fig. IHJl shows the results of UHF calculations 
without correlations, and the top line (offset) in Fig. |SJl presents values of A2 obtained by taking correlations into 
account. Fig. [SJa presents the excitation gap with correlations included. As in the round dot, including correlation 
effects reduces the ground state spin in several cases and hence leads to the violation of Hund's rule. In this case the 
excitation gaps for these N can be larger, yet the ground state spin still changes to the minimal value. However, most 
of the A2 peaks given by the shell structure are still visible. The UHF results of Fig. ^p are qualitatively similar, 
although not identical, to the mean- field (density functional theory) results of Austing et al7 for the Wxjijjy — 1.1 
dot. 

Fig. I^shows the same quantities A2 and A for for dot with large eUipticity u^xl'^y — 2 and with Wat = J^x^y = 0.73. 
The range for r^ over N remains the same, r^ = 1.48. . .0.99. As in Fig. [7^, the spin identification within UHF is 
ambiguous and should be treated with caution. The ambiguity disappears when correlations are included. It is known 
that there are shells in a parabolic potential with LOx/i^y = 2 (see Ref.^^). However the self-consistent potential is 
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not parabolic and we do not observe a shell structure in the dot, with or without correlations. This is evident from 
the gap values plotted in Fig. ^. However, there is some peculiarity at iV = 17,18, and 19 and this may perhaps 
be considered as reminiscent of the shell structure. The uux/^y = 2 case was also included in the mean-field study of 
Austing et alJ and no clear signature of the shell structure was observed. 

VII. CONCLUSIONS 

We have performed calculations for small quantum dots (N < 22) taking into account electron-electron correlations 
in the second order of perturbation theory. We believe that this approach is well justified for r^ < 3 when the 
correlation correction does not exceed a few percent of the total energy. To check the method we have performed 
a comparison with accurate solutions for two and six electrons, in which the method was valid to beyond rg — 3. 
We demonstrate that correlations are especially important for the spin structure of the dot. In some cases even at 
rs = 1 correlations change the spin of the ground state from that found in Hartree-Fock and spin-density-functional 
calculations. In some situations correlations destroy Hund's rule for open electronic shells. Finally, we observe that 
correlations destroy the static spin-density waves observed in Hartree-Fock and spin-density-functional calculations 
and hence restore the spin rotational invariance of states with total spin zero. 

We thank P. G. Silvestrov and A. Micolich for helpful discussions. We are also grateful for the computing resources 
provided by the Australian Partnership for Advanced Computing (APAC) National Facility. 
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